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Summary: Calculations for mean field dynamo models (in both full spheres and 
spherical shells), with both algebraic and dynamic a-quenchings, show qualitative as well 
as quantitative differences and similarities in the dynamical behaviour of these models. We 
summarise and enhance recent results with extra examples. 

Overall, the effect of using a dynamic a appears to be complicated and is affected by the 
region of parameter space examined. 
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1. INTRODUCTION 

In most studies of axisymmetric mean field dynamo models (see for example, 
Tavakol et al., 1995), the nonlinearity is introduced through an algebraic form of a- 
quenching. Such models have produced a number of modes of behaviour in spherical 
and spherical shell dynamo models, including periodic, quasiperiodic and chaotic 
solutions. In addition it has recently been shown that spherical shell models with 
this type of quenching are capable of producing various forms of intermittent type 
behaviour (Tworkowski et al., 1998), which could be of relevance in understanding 
some of the intermediate time scale variability observed in the output of the Sun 
and stars. 

Since such algebraic forms of a-quenching act instantaneously, it is possible that 
this may have a bearing on the occurrence of the more complicated modes of be¬ 
haviour, such as chaos and intermittency, observed in these models. Here using 
recent results as well as new ones, we make a brief comparison between models with 
dynamic and algebraic a-quenching and examine whether features, such as chaotic 
and intermittent-type behaviour, survive as a-quenching is made dynamic. 

2. THE MODEL 

The standard mean field dynamo equation (cf. Krause and Rddler, 1980) is 

= V X (u X B -I- aB — rytV x B), (1) 

at 
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where as usual, the magnitudes of the a and to effects are given by the dynamo 
parameters and C^. 

In the algebraic case we use a purely hydrodynamical a with a = ah, for which 
we take the usual functional form of a-quenching 


Oo cos 6 

^ 1 + B2' 


( 2 ) 


For the dynamical case, following to Zeldovich et al. (1983) and Kleeorin and Ruz- 
maikin (1982) (see also Kleeorin et al., 1995), a can be divided into a hydrodynamic 
and a magnetic part, 

a = ah + a^n, (3) 

where ah is given by Eq. (2) and the magnetic part satisfies an explicitly time 
dependent diffusion type equation with a nonlinear forcing in the form 


dam 
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( 4 ) 


(see Covas et al, 1997a,b; 1998a). This equation differs from that used by Kleeorin 
et al. (1995) in that we take the damping term to be that of the form VaKI^am 
instead of —Om/T, where T is some damping time. Our approach is motivated by 
stability considerations. 

We solved the above equations using spherical polar coordinates. We shall con¬ 
sider both axisymmetric spherical and spherical shell models, where the outer bound¬ 
ary in both cases is denoted by R and in the spherical shell models, the fractional 
radius of the inner boundary of the shell is denoted by r^. We discuss the behaviour 
of the dynamos by calculating the total magnetic energy, E, which is split into two 
parts, E = iffol -b where E^^l and E^^l are respectively the energies of the 
field whose toroidal component is antisymmetric and symmetric about the equator. 
The overall parity P given by P = (Pl®l — E^^l)(E, with P = —1, -bl denoting the 
antisymmetric (dipole-like) and symmetric (quadrupole-like) pure parity solutions 
respectively. 

For B we assume vacuum boundary conditions and for a^, we use a^, = 0 on the 
inner and outer boundary. 


3. RESULTS 

For the the algebraic a-quenching model, we solved equations (1) and (2), whilst 
for the dynamical case we solved equations (1) and (4). Clearly our conclusions 
are subject to the finiteness of the resolution of parameter space that we chose. 
We also point out that the parameters Ca and do not play identical roles in 
these models and therefore the comparison of the behaviours in the two cases is not 
straightforward. 
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3.1 Spherical dynamo models 

Taking = —10"^, which is a typical value that has been used in other dynamo 
studies, for example Tavakol et al. (1995) (and references therein), we found that as 
Ca became large, the algebraic model tended to an antisymmetric (dipolar) parity 
whilst the dynamic case tended towards a symmetric (quadrupolar) one. Overall, no 
chaos was observed in either case although the dynamic case did show evidence of 
more complicated behaviour in its transitions from one parity form to another. We 
also wished to study the allowed forms of behaviour in the supercritical regimes, so 
we chose values which were effectively the highest values numerically allowed by 
our code, and these turned out to be —lO"' and —10® in the dynamic and algebraic 
cases respectively. For the algebraic model, we found chaotic behaviour at large 
value when = —10® was used. 

3.2 Spherical shell dynamo models 

Using C„ = —10"^, we found that for the dynamical a-quenching case, the be¬ 
haviour observed for large Ca values depended on the shell thickness. For thick shells 
(ro = 0.2) we observed symmetric parity solutions whilst for medium (ro = 0.5) and 
thin (ro = 0.7) shells we found antisymmetric asymptotic parity solutions. In par¬ 
ticular, for the medium shell we observed multiple attractors at Ca ~ 15, that is, we 
observed the occurrence of different dynamical solutions possessing different mean 
energies at the same value of Ca but having different initial parities. Also, chaotic 
behaviour was observed in the thin shell at Ca values around 35. However, the alge¬ 
braic a-quenching model did not exhibit any chaotic behaviour nor the presence of 
multiple attractors at this value of The asymptotic behaviour obtained for the 
algebraic case was a symmetric parity solution, an oscillatory mixed parity solution 
which varied between symmetric and antisymmetric values but was periodic and a 
non-oscillatory mixed parity solution for the thick, medium and thin shells respec¬ 
tively. However, when was set to —10®, thus in the supercritical regime for this 
case, the asymptotic behaviour was chaotic for all shell thicknesses. 

Our results appear to indicate that dynamic a-quenching smoothes out chaotic 
behaviour possibly because of the back reaction of a^ via diffusion. 

4. INTERMITTENCY 

Spherical shell dynamo models with algebraic a-quenching have been shown to be 
capable of producing various forms of intermittent-type of behaviour (Tworkowski 
et al., 1998), in which the dynamical modes of behaviour for which the statistics 
taken over different time intervals are different. One example of such behaviour 
is given in Figure 1 (see Tworkowski et al., 1998 for other examples). As can be 
seen, for times around 20 to 60 units, the parity is almost antisymmetric but that 
there are also time intervals which are interrupted by excursions of the parity to 
values well away from —1. We have also searched for occurrences of intermittency in 
the dynamic a-quenching model using the algebraic a-quenching given by equation 
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Fig. 1. Total energy and parity for the a-profile given by (2). Cq = 1.94, Vg = 0.4, F = 0 


(2) with no success, at least to the resolution of our parameter regime. However, 
we did find examples of intermittency using another form of algebraic quenching, 
namely, that due to Kitchatinov (1987). This form was derived in the context of 
A-quenching by rapid rotation and is essentially an interpolation formula having 
the correct asymptotic behaviour. An example of the intermittent behaviour found 
using this form is given in Figure 2. Although these two examples appear similar in 
that they both have intervals where the parity is nearly antisymmetric, this does not 
mean that these two are examples of the same type of intermittency. To resolve this 
issue, one needs to examine the appropriate theoretical framework for the detailed 
properties of the corresponding model such as scaling laws (see, for example, Ashwin 
et al., 1998; Covas et al., 1998b). 

Overall, then, it appears that dynamic a-quenching appears to suppress the ex¬ 
istence of intermittency, at least to the resolution of our parameter search. 

5. CONCLUSION 

Our studies indicate that in the full sphere models, the main similarities are 
the presence of similar modes of parity behaviour and the absence of intermittency, 
whilst the differences are in the details of the transitions between the different modes 
of behaviour and the presence of chaotic behaviour in the algebraic a-quenching case. 
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Time t 


Fig. 2. Result for ro = 0.5 with the dynamic a model and a particular form of an (see 
text). The parameters are Ca = 9.34, C„ = —10'* and F = 0. 


For the spherical shell models, we observe differences which depend on the region 
of the parameter space considered. The dynamic a models appear to produce more 
varied modes of behaviour for C^, = —10"*. However, taking the numerical upper 
bound of in each case appears to indicate that the introduction of dynamic a- 
quenching drastically reduces the likelihood of the occurrence of chaotic behaviour 
and intermittency, which was observed in models with algebraic a-quenching. Our 
dynamic a-quenching models also show multi-attractor regimes with the possibility 
of the final state sensitivity (fragility) with respect to small changes in the initial 
parity. 

With regards to the extra complexity introduced by a time dependent a-quen¬ 
ching, our present results show that the behaviour seen in the solutions is complicated 
and with the outcome depending on the region of parameter space considered, rather 
than a simplistic decrease or increase in complexity. 
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